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The decays K^ s — > 7r°7r°z/P and K^^s ~^ tt^tt^vu involve the weak transi- 
tion from a strange to a non-strange quark. Although they have considerably 
■ smaller branching ratios than those for the corresponding rare processes in- 

volving single pions, some may be more distinguishable experimentally from 
, background processes and thus could provide another probe of the s — > dvv 

\^ \ transition. Using recent knowledge of the CKM matrix elements and measure- 

' ments of related processes, we give improved predictions for their branching 

^ ■ ratios both within and beyond the Standard Model. 
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I. INTRODUCTION 



Historically, rare kaon decays have provided a crucial testing ground in which to study 
flavor- changing neutral current (FCNC) and CP-violating phenomena. Within the context 
of the Standard Model, the observation of such decays gives us further information on 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements and provides independent checks 
on the consistency of our understanding drawn from other measurements, such as those of 
semileptonic decays and B-meson decays. They also serve as a good place to look for physics 
beyond the Standard Model. Even if new physics is first observed outside the kaon sector, 
we will want to know the footprint it leaves in rare kaon decays and specifically on the s ^ d 
weak transition. 

In this paper we focus on the decays K^ s ~^ 7r°7r°z/z/ and K^ s ~^ 7i~^'n-~uu. As noted 
previously , at the quark level these FCNC processes involve the s dvv transition. In 
the Standard Model, this transition is dominated by short-distance contributions involving 
loop diagrams that contain W and Z bosons and heavy quarks 0. With the top quark 
mass much greater than that of the charm quark, the imaginary part of the amplitude for 
this transition arises almost entirely from loop diagrams with top quarks and the resulting 
amplitude is then proportional to the CKM factor ImiVt/Vts)- 

The decays, -k^vv and Kl tt^uu are both governed by the same s dvv transi- 

tion. The latter decay is CP-violating and dominated to high accuracy by the short-distance 
contribution. As noted above, its amplitude is consequently proportional to ImiVtdVts) in 
the Standard Model. For both decays, the hadronic matrix elements of the relevant weak 
current can be related to ones which enter charged current semileptonic decays, whose di- 
rect measurement bypasses any theoretical uncertainty in the hadronic matrix element as 
well. The branching ratios that are predicted in the Standard Model lie between roughly 
10^^^ and 10^^°. With negligible long-distance contributions and little hadronic uncertainty, 
these decays have been pointed to as crucial for precision experimental tests of the Stan- 
dard Model and correspondingly as places to look for new physics if the Standard Model 
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fails Because of this, they are being pursued experimentally in spite of very difficult 

experimental backgrounds. 

The rare decays K^^s ~^ ttttuu that we study in this paper have the same advantage of 
allowing a theoretical clean study of the s — > duu transition, although some of the relevant 
semileptonic matrix elements are not as accurately measured experimentally. As shown 
previously in the Standard Model, they unfortunately have the serious disadvantage 
that their predicted branching ratios are several orders of magnitude smaller than for the 
decays involving a single pion. It is nevertheless interesting to pursue them because they 
involve different combinations of the CP-conserving and CP-violating parts of the s — > duu 
transition, and they thus provide additional handles on both the real and imaginary parts 
of amplitude. Furthermore, from an experimental point of view, some or all of these decays 
may prove to be more susceptible to the extraction of a signal from the background. 

In this paper, we significantly refine the predictions for K ttttuu decays in the Stan- 
dard Model using recent knowledge of the CKM matrix elements and measurements of decay 
rates for related processes. We also examine how large these branching ratios could be for 
physics beyond the Standard Model and find that there are significant, model-independent 
limits from other measurements. The paper is organized as follows: In Section 2, we pro- 
vide the general framework for studying the s — duu transition and the contributions to 
it both within and beyond the Standard Model. Section 3 sets forth its relationship to the 
K -rniuu decays under discussion, the CP properties of the amplitudes involved in those 
decays, and the connection of the relevant hadronic matrix elements to measured semilep- 
tonic decays. Numerical results are given in Section 4, followed by some conclusions in 
Section 5. 



II. GENERAL FRAMEWORK 



The effective Hamiltonian for s 



dvv transitions takes the form 



v/2 27r sin^ 9w 



(1) 



3 



where the short- distance physics is lumped in Wds- In the Standard Model, one-loop contri- 
butions to Wds are dominated by penguin and box diagrams with intermediate charm and 
top quarks: 



wir = K,X{x,) + KdX{x,) , (2) 

where A*^ = V*. Vid, with Vij the appropriate CKM matrix element, and Xi = ffil/Mw'^- 
The QCD corrections to the short-distance contributions X{xi) have been calculated some 
time ago in leading order [^,0 and then in next-to-leading order Since the top-quark 

mass is comparable to the weak scale, these corrections are very small for as can be 

seen explicitly in the values given p],^ for X{xt) when written as X{xt) = ?7tXo(x(), with 
the QCD-uncorrected top quark contribution B 



Xt + 2 3xt-6 

+ 7 7T? ^og{xt) 



(3) 



_Xt - 1 {xt 

and the QCD correction factor rjt = 0.994. On the other hand, these corrections have 
considerable importance for X{xc)- 

The quantity X{xt) is roughly three orders of magnitude larger than X{xc), and since 
ImA^^ = — ImA*^, the top contribution completely dominates in the imaginary part of WfJ^. 
However, ReA*^ << ReA^^, allowing the charm contribution, although still smaller in mag- 
nitude than that from top, to be roughly comparable and to interfere constructively in the 
real part of W^J^ . 

As illustrative examples of physics that lies beyond the Standard Model, we consider two 
very different possibilities: 

• Effective Flavor-Changing Neutral Current (FCNC) interaction. Such an interaction. 



as formulated by Nir and Silverman p!0|JTl[| , takes the form of an extra term in the 
effective Lagrangian of the form: 

4 cos Ow 
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When combined with the couphng of the Z boson to neutrino-antineutrino pairs, one 
finds that 



TT 
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as the new piece of Wds in the effective Hamiltonian that corresponds to the basic 
process, s dvv. 

Upper bounds for Uds have been determined by other processes involving K mesons 
and were summarized recently to be 



|Re(t/d.)| < 10-^ (6) 

\Uds\< 3 X 10-^ (7) 

|Re([/rf,)Im(f/d,)| < 1.3x 10-^ (8) 

|Im(f/d,)| < 10-5. (9) 

The bound on \Uds\ arises from the decay vr+z/z/, whose width is proportional 

to It can be improved by using the most recent measurement [|13[, of the 



branching ratio, BR(i^' 10 This value is consistent with what 

is expected in the Standard Model and corresponds to \Wds\ = O.OSigl^ x 10"^. If we 
were to assume that the total branching ratio were due to new physics arising from Uds, 
then the bound on \Uds\ would be reduced from that in Eq. (7) to \Uds\ < 1.6 x 10^^. 

Supersymmetry. A dominant supersymmetric effect arises from penguin diagrams 
involving charged-Higgs plus top-quark intermediate states or squark and chargino 
intermediate states. These give additional pieces to the effective Hamiltonian of the 



form |]T1 



where tan/? is the ratio of the two Higgs vacuum expectation values and XtH = 
rrit^ /Mh±'^- The quantity H{x) is given by 



The parameter Xt can be bounded by similar considerations to those that were used 
for Uds- The observed branching ratios for the decays Kl fi^fi~ and tx^vv 



have been used to set the hmits 1113 



iReAtI < 0.21 , (12) 
\\t\ < 0.35 . (13) 

The most recent branching ratio for ix^vv could be used to revise the last limit 

to \\t\ < 0.16 . 

As we will see shortly, the limitations imposed by experiments on the parameters of both 
these examples of physics beyond the Standard Model lead to similar restrictions on how 
large the branching ratios can be for the processes we are studying. 



III. K -K-KvD DECAYS 

When the effective four-fermion operator relevant for the decay we are considering is 
sandwiched between the initial and final states, it factorizes into a product of matrix elements 
of the hadronic current and the leptonic current. We will use isotopic spin to relate the 
hadronic matrix elements relevant to Kl s ~^ ttttz/z/ to those for -r-kc^v, where the 

corresponding branching ratios (and hence squares of matrix elements) have been measured. 

We consider first the process Kl — * 7r°7r°z/z/. The vr^vr" pair forms a CP-even state, and 
has total isospin, J = 0. The vv pair, created by a virtual is CP even as well. Since 
there must be one unit of orbital angular momentum to allow the total angular momentum 
of the final state to be that of the initial Ki, namely zero, the final state is CP-odd. The 
overall decay process is then CP-conserving for the major (CP-odd) piece of the Kl, and 
the resulting amplitude is proportional to the real part of Wds- The opposite result holds for 
the piece of the Kl that is proportional to e and CP-even; the corresponding decay process 



is CP- violating and the amplitude is proportional to the imaginary part of Wds- It is thus 
suppressed on two counts and is completely negligible. 
Using the relationship, 

(ttV" I (st/)v_A I K") = (vrV" I (sm)v-a I K+) , (14) 

we find that 



BR{K, Mu) = ^Y.^T^'S ,, ^mK^ - vr-TT^e+z.) , (15) 

27i^sm*9w\Vus\ tr^ 



where the factor of 3 accounts for the three species of neutrinos. 

By relating the desired branching ratio to a measured one, we have avoided having either 
to do a calculation of the hadronic matrix elements or to perform a detailed analysis in terms 
of invariant amplitudes, as was done in previous analyses Of course, the final results 

for the branching ratio must be consistent, since both approaches agree with the available 
data on charged-current semileptonic decays, and in particular those for the decay rate for 
the process 'K°'K°e^v. For the purposes of this paper of discussing the absolute and 

relative size of the various branching ratios within and beyond the Standard Model, it is 
considerably easier to formulate the results directly in terms of relationships to branching 
ratios for measured semileptonic decays. 

A similar formula can be obtained for BR{Ks —>■ n^n^jyu), but with KeWds replaced 
by ImVFrfs and tkj^ replaced by rxg- Since the Ks has a much shorter lifetime and the 
major part of the Ks corresponds to a transition that is CP-violating, this branching ratio 
is orders of magnitude smaller than that for 'k°tx°vi'. Although the part of the Ks 

state proportional to e corresponds a CP-conserving transition, the smallness of e still gives 
rise to a net decay amplitude that is much smaller than that for the CP-even part of the 
Ks- 

The analysis for the decay Kl ti^ti^uu can be carried out analogously. It is convenient 
to break it up into the cases where the tt+tt" pair in the final state has total isospin zero 
and one, since there is no interference between them in the decay rate. For the isospin zero 



case, the argument about the CP properties of the final state is the same as given before, 
and we simply have a factor of two in the rate for the tt^tt^ final state compared to that for 
the 7i°7[" final state discussed above: 

BRiKL - (7r+7r-)i=oi^^) = - nWu). (16) 

vr^sm Ow\VusrTK+ 

This situation is slightly more complicated for the case where the tt+tt^ pair has isospin, 
7 = 1. The TTTT pair is still CP-even, but it must be in a p-wave. There are two possible ways 
in which the total angular momentum of the final state can be zero, which correspond to the 
relative orbital angular momentum of the vrvr and uu pairs being zero or one. As is expected 
when there is such limited phase space, the latter amplitude is strongly suppressed by 
centrifugal barrier effects compared to the former [0] . So we are left with a single amplitude 
where the relative orbital angular momentum is zero. The CP of the final state is even and 
the transition involving the major part of the Kl is CP-violating. Using the relationship, 

V2 ((7r+7r-)i=i I (sd)v-A I K") = ((vr-7r°)i=i | (sn)v-A I K") , (17) 

we find that 

BR{Kl (7r+7r-)i=iz/j^) = , , \, .. BRiKL - nVe^u) . (18) 

47r^ sm fc'iyl Ksl 

The total branching ratio for ti^ti^uu is then found by simply adding the two results 

above: BR{Kl n+n-uu) = BR{Kl {n+n-)i=oJ'u) + BR{Kl (7r+7r-)i=iZ/z/). 

The corresponding formula for Ks ti^tt^uu can again be obtained by the interchange 
of ReWsd and ImW^d and multiplication of the right-hand-side by tksI'TKl- 

IV. NUMERICAL CALCULATION 

To obtain numerical predictions we have used a set of parameters taken from the Review 
of Particle Physics |T5[, including the fine-structure constant at the weak scale, a = 1/129; 



Mw = 80.3 GeV; sin^ 9w = 0.23; and the measured semileptonic branching ratios needed in 
Eqs. (15), (16) and (18). We correspondingly find that: 
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BR{Kl ■K°Ti°vv) = [(3.1 ± 0.6) X 10"^] \ReW, 
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ds I 

BR{Kl {n+n-)i=oi^i^) = [(6.2 ± 1.2) x 10"^] \ReWds\^ , 

BR{Kl (7r+7r-)i=iz/z/) = [(0.93 ± 0.05) x 10"^] \lraWds\'^ , (19) 

where the error bars come from those of the experimental measurements of the relevant 
semileptonic branching ratios. We have not taken account of radiative corrections or isotopic- 
spin violating differences in form factors and phase space, as has been done for the case of 
the decays involving a single pion |jl6|, given the (larger) uncertainties in other parts of 



input at this stage of the analysis of these decays. Formulas similar to Eq. (19) hold for the 
related decays of the Ks to the same final states. 

For the specific calculation of Wds in the Standard Model we need the values of X{xt) 
and of X{xc) in next-to-leading order [P,|l, and that of the CKM matrix elements flj 



ReVtd = 0.0076 ± 0.0015 , 

ImVtd = 0.0031 ± 0.0008 , (20) 

and Vts = —Vcb = —0.040 ± 0.002, aside from the well-known matrix elements connecting 
the first and second generations. Using this and with rrit = 166 ± 5 GeV, we find that 

= [(-6.7 ± 1.0) + i{1.9 ± 0.5)] X 10"^ , (21) 

and the branching ratios for the various processes shown in Table I. 

In Kl decays the contribution of the isospin one tt+tt^ final state is negligible, since it 
is already suppressed compared to that with isospin zero from Eq. (19) and the magnitude 
of the real part of W^J^^ is considerably greater than that of the imaginary part. Thus the 
ratio between tt+tt^ and n"n° rates is very close to the factor of two characteristic of isospin 
zero. 

Our results are given in Table I, and both Kl branching ratios lie between 10~^^ and 
10^^^. Our predictions in the Standard Model for BR{Ki 7r°7r°z/z/) are consistent with the 
previous calculation [0 of 1 - 3 x lO^^^ and those [0,0 of 1.1 - 5 x lO^^^ and 2 - 5 x lO^^^ 



Processes 




Br(x 10-^3) 


SM 


New Physics (maximal) 




1.4 ±0.4 
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Ks TT^TT^VV 


(1.9 ±0.8) X 10"^ 


8 X 10-^ 


Kl — ^ 'k'^'k^uv 


2.8 ±0.8 


20 




(11 ±2) X 10""^ 


2 X 10-2 



TABLE I. Branching ratios of various Kl,s TrTrvv modes within and beyond the Standard 



Model. 

for BR(Kl — > TT+TT-i/p), but the allowed range is now considerably restricted. The Ks 
branching ratios are in the 10"^^ range in the Standard Model. The decay Ks t:^t:~'ui> 
gets important contributions from both 7 = and 7 = 1 tttt final states since the suppression 
of the 7 = 1 final state is compensated by the ratio of {ReWds/^'^WdsY in the Standard 
Model. 

For the representative examples of physics beyond the Standard Model, we also show 
in Table I values for the branching ratios that correspond to the maximal values one could 
obtain consistent with the bounds in Eqs. (6)- (9) and (12) - (13), respectively. These 
maximal values are similar in both cases and arise when the parameters of the new physics 
are chosen to maximize W^s^ consistent with the constraints coming from known K physics. 
Among these constraints is the recent branching ratio measurement for K^ n'^uu, which 
is equally sensitive to both the real and imaginary parts of the total Wds with minimal 
theoretical assumptions. Hence, similar maximal values of W^^^ are obtained in any model 
of new physics. Note that when taking ratios to the Standard Model branching ratios, a 
much bigger factor is possible when the new physics enters the imaginary part of W^^, and 
is CP- violating, since the imaginary part of W^J^ is considerably smaller than the real part. 
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V. SUMMARY 



We have used recent information on the CKM matrix to narrow the range of the predicted 
branching ratios for K^ s ~^ 'k°'k°i'D and K^^s ~^ 7i~^7i~i'iy decays in the Standard Model. 
These branching ratios are in the neighborhood of 1 — 4 x 10~^^ for the Kl decays, which 
make them possibly observable at the few event level in the next round of experiments that 
are setting out to see the CP-violating decays with a single 7r° in the final state. These 
branching ratios could be larger by up to about an order of magnitude in theories that go 
beyond the Standard Model. 

The branching ratios for the Ks decays are in the neighborhood of 10~^^ to 10~^^ in the 
Standard Model and seem unlikely to ever be observed. Here new physics could boost the 
branching ratios by more than an order of magnitude, although even then the maximum 
branching ratio of around 10~^^ is still beyond the limits of observation. For both Ks 
and Kl decays, increased experimental accuracy in the measurement of the branching ratio 
for K'^ — > TT'^vP, assuming it remains consistent with the Standard Model, will put more 
stringent restrictions on non-Standard-Model physics in the s dvP transition and limit 
the deviations from the Standard Model that can be observed in the decays under discussion 
here as well. 
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